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Abstract 

We discuss recent calculations of the survival probability of the large rapidity gaps in 
exclusive processes of the type pp — > p + A + p at high energies. Absorptive or screening 
effects are important, and one consequence is that the total cross section at the LHC is 
predicted to be only about 90 mb. 

At the LHC, the observation of an exclusive process of the type pp — > p + A + p, where a 
produced new heavy object A is separated from the outgoing protons by large rapidity gaps 
(LRG), will provide very good experimental conditions to study the properties of object A 
[1, 2, 3]. The process is sketched in Fig. 1. The case of A = H —>■ bb is particularly interesting. 
Thecross is usually written in the form 



(S 2 ) 
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where B/2 is the t-slope of the proton- Pomeron vertex, and the constant is known in 
terms of the A — ► gg decay width. The amplitude-squared factor, |...| 2 , can be calculated 



l Based on a talk by A.D. Martin at the CERN - DESY Workshop "HERA and the LHC" ,26 - 30 May 2008, 
CERN. 




Figure 1: The mechanism for the exclusive process pp — > p + A + p, with the eikonal and 
enhanced survival factors shown symbolically. 

in perturbative QCD, since the dominant contribution to the integral comes from the region 
Aq CD <C Ql <C M A , for the large values of M\ of interest. The probability amplitudes, f g , to 
find the appropriate pairs of t-channel gluons (xi,x[) and (x 2 , x' 2 ) of Fig. 1, are given by skewed 
unintegrated gluon densities at a hard scale \i ~ Ma/2. To evaluate the cross section of such 
an exclusive processes it is important to know the probability, (S 2 ), that the LRG survive and 
will not be filled by secondaries from eikonal and enhanced rescattering effects. The main effect 
comes from the rescattering of soft partons, since they have the largest absorptive cross sec- 
tions. Therefore, we need a realistic model to describe soft interactions at the LHC energy, and 
to predict the total cross section at LHC. The model must account for (i) elastic rescattering 
(with two protons in intermediate state), (ii) the probability of the low-mass proton excitations 
(with an intermediate proton replaced by the N(1400), N(1700), etc. resonances), and (hi) the 
screening corrections due to high-mass proton dissociation. 

The effect of elastic rescattering may be evaluated in a model independent way once the 
elastic pp-amplitude is known. The effect of the low-mass dissociation is usually calculated in the 
framework of the Good- Walker formalism[4], that is, by introducing diffractive eigenstates, 0« 
with % = 1, ..n, which only undergo 'elastic' scattering. The resulting n-channel eikonal fljfc(s, b) 
depends on the energy and the impact parameter of the pp interaction. The parameters of the 
model are chosen to reproduce the available (fixed-target and CERN-ISR) data on the cross 
section of low-mass diffractive dissociation. Usually either a two- or three-channel eikonal is 
used. Finally, high- mass dissociation is described in terms of Reggeon diagram technique [5]. A 
symbolic representation of these soft scattering effects is shown in Fig. 2. The latest calculations 
along these lines are described in Refs. [6, 7]. In Ref. [6] the authors account only for the triple- 
Pomeron vertex, and, moreover, sum up only the specific subset 2 of multi-Pomeron diagrams 
that were considered in Ref. [8], which is called the MPSI approximation. In Ref. [7] all 

2 For example, the third, but not the second, term on the right-hand side of the expression for f2jfc/2 in Fig. 2 
is included; neither are multi-Pomeron terms, like the last term, included. 
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Figure 2: The multi-channel eikonal form of the amplitude, where i, k are diffractive (Good- 
Walker) eigenstates. Low-mass proton dissociation is included by the differences of the Pomeron 
couplings to one or another Good- Walker state (i) in the first diagram, while the remaining 
(multi- Pomeron) diagrams on the right-hand side of the expression for 1^/2 include the high- 
mass dissociation. 

possible multi-Pomeron vertices were included under a reasonable assumption about the form 
of the n —>■ m multi-Pomeron vertices, g^. The assumption corresponds to the hypothesis that 
the screening of the s-channel parton c during the evolution is given by the usual absorption 
factor exp(—Q ic (b) — Q c k(b)), where VL ic {b) (ft c k(b)) is the value of the opacity of the beam 
(target) proton at impact parameter b with respect to the parton c. 

Since the absorptive corrections increase with energy, the cross section grows more slowly 
than the simple power (er oc s A ) parametrisation [9]. In spite of the fact that the models of [6] 
and [7] are quite different to each other, after the parameters are fixed to describe the data on 
the total, elastic and single dissociation cross sections (cr to t, da e \/dt and dasn/dM 2 ) within the 
CERN-ISR - Tevatron energy range, the latest versions of the Tel-Aviv and Durham models 
predict almost the same total cross section at the LHC, namely a tot ~ 90 mb. Correspondingly, 
both models predict practically the same gap survival probability (£f ik ) ~ 0.02 with respect to 
the eikonal (including the elastic and low-mass proton excitation) rescattering, for the exclusive 
production of a Higgs boson. 

A more delicate problem is the absorptive correction to exclusive cross sections caused by 
the so-called enhanced diagrams, that is by the interaction with the intermediate partons, see 
Fig. 1. This rescattering violates 'soft-hard' factorisation, since the probability of such an 
interaction depends both on the transverse momentum and on the impact parameter of the 
intermediate parton. 

The contribution of the first enhanced diagram was evaluated in[10] in the framework of 
the perturbative QCD. It turns out to be quite large. On the other hand, such an effect is not 
seen experimentally. The absorptive correction due to enhanced screening must increase with 
energy. This was not observed in the present data (see [11] for a more detailed discussion). 

Several possible reasons are given below. 

(a) We have to sum up the series of the multi-loop Pomeron diagrams. The higher-loop 
contributions partly compensate the correction caused by the first-loop graph. 
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(b) There should be a "threshold", since Pomeron vertices must be separated by a non- 
zero rapidity interval [12]. That is, at present energies, the kinematical space available for 
the position of a multi-Pomeron vertex in an enhanced diagram is small, and the enhanced 
contribution is much less than that obtained in leading logarithmic (LL) approximation. 

(c) The factor S^ ik already absorbs almost all the contribution from the center of the disk. 
The parton only survives eikonal rescattering on the periphery, that is at large b. On the other 
hand, on the periphery the parton density is rather small, and the probability of enhanced 
absorption is not large. This fact can be seen in Ref. [13]. There, the momentum, Q s , below 
which we may approach saturation, was extracted from HERA data in the framework of the 
dipole model. Already at b = 0.6 fm the value of Q 2 S < 0.3 GeV 2 for x < 10~ 6 . See also [14] 
where the value of Q s was evaluated using LO DGLAP evolution. 

Point (c) is relevant to the calculation of S^ nh described in [6]. First, note that the b 
dependence of the beginning of 'saturation caused by enhanced graphs' is not accounted for in 
the MPSI approximation used in [6]. In this model, we have the same two-particle irreducible 
amplitude (which sums up the enhanced diagrams) at any value of b. Therefore, the enhanced 
screening effect does not depend on the initial parton density at a particular impact parameter 
point b. For this reason the suppression due to enhanced screening corrections (S^ nh ) = 0.063 
claimed in [6] is much too strong 3 . 

The survival factor (S^ nh ) has also been calculated in the new version of the Durham model 
[16]. The model includes 3 components of the Pomeron, with the different transverse momenta 
k t of the partons in each Pomeron component, in order to mimic BFKL diffusion in In k t . In 
this way we obtain a more realistic estimate of the 'enhanced screening' in exclusive diffractive 
Higgs boson production at the LHC. The model predicts (S^ nh ) ~ 1/3. However the CDF data 
on exclusive 77 and Xc production indicate that this suppression is not so strong. 

Note, that comparing the values of the survival factors in this way is too simplistic. The 
problem is that, with enhanced screening on intermediate partons, we no longer have exact 
factorisation between the hard and soft parts of the process. Thus, before computing the effect 
of soft absorption we must fix what is included in the bare exclusive amplitude calculated in 
terms of perturbative QCD. 

The first observation is that the bare amplitude is calculated as a convolution of two gen- 
eralised (skewed) gluon distributions with the hard subprocess matrix element, see (1). These 
gluon distributions are determined from integrated gluon distributions of a global parton anal- 
ysis of mainly deep inelastic scattering data. Now, the phenomenological integrated parton 

3 Moreover, since the irreducible amplitude approaches saturation at some fixed energy (rapidity), indepen- 
dent of the value of b, the approximation gives <7tot(s — * 00) — » constant. On the other hand, a theory with 
an asymptotically constant cross section can only be self-consistent in the so-called 'weak coupling' regime for 
which the triple-Pomeron vertex vanishes for zero momentum transfer [15]. The vertex used in [6] does not 
vanish. This indicates that the MPSI approximation cannot be used at asymptotically high energies, and the 
region of its validity must be studied in more detail. 
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distributions already include the interactions of the intermediate partons with the parent pro- 
ton. Thus calculations of S en h should keep only contributions which embrace the hard matrix 
element of the type shown in Fig. 1. 

The second observation is that the phenomenologically determined generalised gluon distri- 
butions, f g , are usually taken at p t = and then the observed "total" cross section is calculated 
by integrating over p t of the recoil protons assuming the an exponential behaviour e~ Bp *; that 
is 



However, the total soft absorptive effect changes the p t distribution in comparison to that for the 
bare cross section determined from perturbative QCD. Thus the additional factor introduced 
by the soft interactions is not just the gap survival S 2 , but rather the factor S 2 /B 2 [17], which 
strictly speaking has the form S 2 (p 2 ) 2 . 

In order to compare determinations of the suppression due to absorptive effects we should 
compare only the values of the complete cross section for pp — > p+A+p. However a comparison 
is usually made by reducing the cross section to a factorized form. If this is done, as in (1), then 
the Durham predictions for the survival factor to eikonal and enhanced screening of the exclusive 
production of a 120 GeV Higgs at the LHC are (S 2 ) = 0.008, 0.017,0.030 where enhanced 
sreening is only permitted outside a threshold rapidity gap Ay = 0, 1.5, 2.3 respectively. The 
values correspond to B = 4 GeV~ 2 . 

Let us discuss the survival factors claimed by Frankfurt et al. [18]. They use another 
approach. Within the eikonal formalism, they account for elastic rescattering only. The pos- 
sibility of proton diffractive excitation is included in terms of parton-parton correlations, for 
both low- and high-mass dissociation. At a qualitative level, it is possible to consider all the 
effects discussed above in terms of such a language. On the other hand, to the best of our 
knowledge, they did not describe the available data on er tot , da c \/dt, M 2 da^ / dM 2 . Also, the 
energy (i.e. 1/x) dependence of the parton densities was evaluated using simple LO DGLAP 
evolution. This is grossly inadequate for the low values of x sampled, x ~ 10~ 5 . Thus, it is 
difficult to judge the accuracy of their numerical predictions. Moreover, part of the Sudakov- 
like suppression, which above was calculated using perturbative QCD, is here treated as parton 
correlations and included in the value of S 2 nh . 4 Therefore, one cannot compare literally the 
predictions for the gap survival factors S* 2 = (S 2 ik (b) S 2 nh (b)) given by [18] and by the Durham, 
Tel- Aviv and Petrov et al. [19] models 5 . The only possibility is to compare the predictions for 
the final exclusive cross section. Unfortunately, such a prediction is not available in [18]. 

4 In general, one may include the absence of QCD radiation in the large rapidity gap in the "soft" survival 
factors, but to make comparisons we must define precisely in which part of the calculation each effect is 
included. Note also that in [18] the DL expression for Sudakov T-factor is used, which grossly overestimates 
the suppression. 

5 The last group calculated S 2 within their own eikonal model and fitted the parameters in a Regge-type 
expression for f g to describe HERA data. The final prediction is again rather close to that by the Durham 
group. 




(2) 
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Next, we comment on another recent calculation [20] along the lines of eq. (1). They claim 
very large uncertainties in the predictions arising mainly from the freedom in the choice of 
limits of integration in the Sudakov form factor which is embedded in f g . However, this is 
not the case. In fact, the Sudakov factors have been calculated to single log accuracy. The 
collinear single logarithms are summed up using the DGLAP equation. To account for the 
'soft' logarithms (corresponding to the emission of low energy gluons) the one- loop virtual 
correction to the gg — > A vertex was calculated explicitly, and then the scale \x = 0.62 M A 
was chosen so that double log expression for the Sudakov form factor reproduces the result of 
the explicit calculation. Similarly, the lower limit k 2 = Q\ was verified to give the one-loop 
result. It is sufficient to calculate just the one-loop correction since it is known that the effect 
of 'soft' gluon emission exponentiates. Thus double log expression, with \x = 0.62 Ma, gives 
the Sudakov factor to single log accuracy. Also the form used for / 9 's in Ref. [20] contradicts 
the known leading log(l/x) asymptotic behaviour. 

Finally, we discuss a very recent calculation [21] based on the dipole approach. A new 
development is that instead of using a multi-channel eikonal with a fixed number of diffractive 
eigenstates, the authors consider an explicit wave function of a fast hadron (proton, pion) 
and have a continuous integration over the size of the quark-quark dipoles. In this model the 
incoming hadron wave function is approximated by a simple Gaussian. The parameters are 
fitted so as to describe the data on <7 tot , <r e i and F 2 at low x. A shortcoming is that high- 
mass dissociation is calculated separately. Its contribution is not included in the proton dipole 
opacity fl(r,b), for which a simplified asymptotic solution of the BFKL equation was used. 
Moreover, to calculate the gap survival probability, S 2 (b), the b dependence is considered, but 
the dependence of the "hard subprocess" cross section on the dipole size was not accounted for. 
That is, again, the correlation between the saturation momentum Q s and b is lost. Nevertheless, 
the model confirms the observation that the energy dependence of S 2 is not too steep; S 2 at 
the LHC for central exclusive production is only reduced by a factor of about 2.5 to that at the 
Tevatron. Thus, Tevatron data serve as a reliable probe of the theoretical model predictions of 
these production rates. 

In summary, we have briefly discussed various recent calculations of the exclusive process 
pp — > p + A + p at high energy. The value of the cross section when A = (H — > bb) is important 
for the feasibility of using tagged protons to study the Higgs sector via this process at the 
LHC. We have paid special attention to the survival factors of the large rapidity gaps. We see 
no reason to doubt the claimed value, or accuracy, of the existing predictions of the Durham 
model. Recall that these predictions have been checked in many places by comparing with 
the available experimental data on exclusive 77 and high Et dijet production at the Tevatron 
and on exclusive diffractive J/ip production at HERA (see [22, 23] for more details). Since all 
the factors, which enter the calculations, depend rather weakly (logarithmically) on the initial 
energy, there is no reason to expect that the model, which describes the data at the Tevatron 
energy, will be too far from reality at the LHC. 
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